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ABSTRACT 

Despite intense scrutiny, the progenitor system(s) that gives rise to Type 
la supernovae remains unknown. The favored theory invokes a carbon-oxygen 
white dwarf accreting hydrogen-rich material from a close companion until a 
thermonuclear runaway ensues that incinerates the white dwarf. However, sim- 
ulations resulting from this single-degenerate, binary channel demand the pres- 
ence of low- velocity Ha emission in spectra taken during the late nebular phase, 
since a portion of the companion's envelope becomes entrained in the ejecta. 
This hydrogen has never been detected, but has only rarely been sought. Here 
we present results from a campaign to obtain deep, nebular-phase spectroscopy 
of nearby Type la supernovae, and include multi-epoch observations of two 
events: SN 2005am (slightly subluminous) and SN 2005cf (normally bright). 
No Ha emission is detected in the spectra of either object. An upper limit 
of 0.01 of solar abundance material in the ejecta is established from the 



models of iMattila et al.l (120051 ) whic h, w hen coupled with t he mass-stripping 



simulations of Marietta et al.l (l2000l ) and iMeng et al.l (120071 ) . effectively rules 
out progenitor systems for these supernovae with secondaries close enough to 
the white dwarf to be experiencing Roche lobe overflow at the time of explosion. 
Alternative explanations for the absence of Ha emission, along with suggestions 
for future investigations necessary to confidently exclude them as possibilities, 
are critically evaluated. 

Subject headings: binaries: symbiotic — circumstellar matter — supernovae: 
general — supernovae: individual (SN 2005am, SN 2005cf ) — white dwarfs 



^Some of the data presented herein were obtained at the W.M. Keck Observatory, which is operated 
as a scientific partnership among the Cahfornia Institute of Technology, the University of California and 
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Introduction 



Ever since the currently favored single -degenerate, binary ch annel was proposed as the 
progenitor system for Type la supernovae (IWhelan &: Ibenlll973l ). all models of the impact 
of the exploded white dwarf (WD) on the secondary star have indicated that significant 



amounts of solar-abundance materia 



trained in th e ejecta (IWheeler et al 



stripped from the secondary ' s envelope, become en- 



I975I: IPrvxell fc Amettl 



evidence for this material, however, still eludes us (IMattila et al.ll2005l ). and serves as one 



1981 



Chugail Il98d : Ihivne et all Il992l : [Marietta etaP I2OO0I: iMeng et all l2007h . Observational 



Taam fc Frvxelll Il984 : 



Branch et al. 


1995; 


Livio 


2001) 



The most detailed theoretical investigation of the expected amount and distribu- 
tion of strippe d mate rial within a young Type la supernova (SN la) remnant is that of 
Marietta et al.l (120001 ). who studied the problem with two-dimensional numerical simula- 
tions. Four basic progenitor systems were investigated, including three with secondaries 
(a main-sequence, subgiant, or red giant star) close enough for mass-transfer to occur 
through Roche lobe overflow (RLOF), and one containing a secondary (a red giant) do- 
nating material through a strong stellar wind — the symbiotic case. All secondaries were 
given masses of ~ 1 Mq at the time of the explosion and were placed either just within 
(the RLOF cases), or just beyond (the symbiotic case), the limiting distance from the WD 
within which RLOF can occur (i.e. , a/R = 3, whe re a is the orbital separation in units of 



the secondary star's radius, R; see lEggletonlll983l ). Three additional systems, in which a 
main-sequence secondary was placed too far away to experience RLOF (thus rendering it 
an unlikely progenitor system) were also included to establish the scaling between orbital 
separation and amount of stripped material. 



The numerical results of the [Marietta et al.l (12000f) study confirmed predictions from 



earlier analytic work (e.g.. IWheeler et al.lll975l : IChugailll986l ) that substantial material is 
indeed stripped, with the amount ranging from a minimum of 0.15 Mq for a close (a/ R = 3) 
main-sequence secondary, up to nearly the entire envelope (~ 0.5 Mq) for a similarly placed 



red gi ant. Increasing the orbital separation beyond the RLOF limit in the [Marietta et al. 



(120001 ) simulations resulted in a dramatic decrease in the amount of stripped material: 
For a/R = 12, a main-sequence secondary loses only 0.0018 Mq. However, since such 



the National Aeronautics and Space Administration. The Observatory was made possible by the generous 
financial support of the W.M. Keck Foundation. Additional observations were obtained at the Gemini 
Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., un- 
der a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science 
Foundation (United States), the Particle Physics and Astronomy Research Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the Australian Research Council (Australia), 
CNPq (Brazil) and CONICET (Argentina). 
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systems lack an efficient mechanism for mass transfer, they are not considered to be viable 
SN la progenitors. A red giant secondary placed at a similarly distant location, however, 
could potentially donate mass through a strong stellar wind, making the symbio tic case a 



possibility if large orbital separations are required (e.g.. iMunari fc Renzinilll992l ). 



As discussed by lMeng et al.l (l2007l ). one shortcoming of the lMarietta et al.l (120001 ) study 
is their use of standard solar-model stars for the companion, rather than companions whose 
structures have been appropriately modified due to having evolved in a binary system (e.g., 
Eggleton I973I). T o investigate the effect this simplification might have had on the results, 
Meng et al.l (120071 ) use an analytic model to estimate the amount of mass expected to be 
stripped from a variety of evolved second aries. (Their a.nalyti c approach was first tested 
using the unevolved secondaries used by [Marietta et al.l I2OOOI and was demonstrated to 
approximate the results obtained numerically.) The result is that the quantity of material 
expected to be stripped from evolved secondari es is conside r ably l ower than that predicted 
for standard solar-model companions. In fact, iMeng et al.l (120071 ) find that the minimum 
value for systems experiencing RLOF is diminished from 0.15 Mq to only 0.035 Mq. The 
reduction arises primarily from the pre-explosion mass-loss producing a more compact 
companion s tar w' hose material is more difficult to strip than it is in the unevolved case. 
Meng et al.l (120071 ) stress, however, that their new values are really lower bounds on the 
amount of stripped material, since their analytic approach does not consider the thermal 
energy imparted by the ejecta to the companion envelope, which likely serves to heat 
and vaporize a portion of it and thereby increase the amount of stripped material (e.g., 
Fryxell fc Arnettlll98ll : lMattila et al.ll2005l ). Thus, 0.035 Mq serves as a conservative lower 
bound on the expected amount of stripped material resulting from their models. 

The typical velocity of the stripped material is found in all studies to be far slo wer than 
the ~ 10,000 km s~^ veloc i ty that characteri zes the bulk of the iron-rich ejecta (jChugai 



19861 : Marietta et al.l I2OOOI : iMeng et all 120071 ). This has the effect of placing it almost 
entirely in the central region of the supernova remna nt, with the majority of it predicted 
to be moving with a velocity of under 1, 000 km s~^ (iMarietta et al.ll2000l ). The stripped 
material is largely confined to the downstream region behind the companion star, where it 
contaminates a solid angle th at ranges from 66° fo r the main-sequence companion to 115° 
for the red-giant companion (IMarietta et al.ll2000l ). 



The low velocity of the stripped material renders it undetectable when the faster- 
moving, iron-rich eje cta are optical l y thic k. However, detailed radiative transfer calcu- 
lations performed by iMattila et al.l (120051 ) predict that it should become visible at late 
times, when the outer ejecta have thinned out and become transparent enough to reveal 
the slower-moving gas in the central regions. The most prominent expected spectral sig- 
nature of the compa nion star's stripped material is narrow Ha emission in nebular spectra 
(IMattila et al.l 120051 ) . taken more than ~ 250 days after maximum light. The Ha emis- 
sion should be present within ±1,000 km s~^ of Aq = 6563 A but, due to the expected 
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asymmetry in the distribution of the solar-abundance material, could present an Ha profile 
ranging from a very narrow spike to a broader emission line in the observed spectrum. 
Obtaining spectra with high enough resolution to resolve fairly narrow lines is thus a useful 
component of a targeted search for this Ha. 

To date, only a few nebular SN la spectra have been obtained and Ha has never been 
detected, although the majority of the spectra lack the spectral resolution and signal-to- 
noise ratio to place intere sting constraints on the companion. The best limits, by far, come 
from the recent study by iMattila et al.l (120051 ) , which constrains hydrogen-rich material in 
SN 2001el to be ^ 0.03 Mq from a low-resolution (~ 700 km s~^) spectrum obtained 398 
days past maximum light. 

In an effort to expand and improve on earlier work, both in terms of the number of 
objects studied as well as the resolution, sensitivity, and temporal coverage of the spectra, 
we have initiated a program to obtain deep, moderate resolution (< 150 km s~^, or ~ 
3 A, at Ha), late-time spectra of SNe la at multiple epochs using the Keck and Gemini 
telescopes. The first p hase of this pro ject has garnered data on two objects: SN 2005am, 
a slightly sublurninous (ILi et al.ll2006l ) event and SN 2005cf, an SN la of normal brightness 
(jPastorello et al.ll2007l ). We present and analyze our observations in §[2], discuss the results 
in § 131 and conclude in § HJ 



2. Observations and Analysis 

We obtained a total of five deep nebular-phase spectra: Two for SN 2005am and three 
for SN 2005cf . Details of the observations are given in Table [U Following initial pro cessing 
of th e frames]^ we extracted all one- dimensional sky-subtracted spectra optimally (Home 



1986 ) in the usual manner using the apall task within IRAfJ^I Each spectrum was then 



wavelength and fiux calibr ated, and corrected fo r continuum atmospheric extinction and 



telluric absorption bands (IMatheson et al.ll2000l . and references therein). Careful exami- 



nation of the two-dimensional spectra reveals only faint night-sky emission at the nominal 
location of Ha, and so we deem the potential for contamination at these wavelengths due 
to improper background ( "sky" ) removal to be negligible. 

The final spectra are displayed in Figure [H To search for Ha emission, or to place 
limits on the amount of solar- abundance material present in the inner ejecta, we subjected 



^For Gemini observations, the frames were processed using the tasks gprepare, gsreduce, gsflat, and 
gmosaic in the Gemini IRAF package. 

•^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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each spectrum to the following analysis procedure. 

First, we attempted to place each spectrum on an absolute flux scale. Adjustments to 
the original flux levels are necessary since most observations were made with a rather narrow 
slit width compared with the seeing (see Tabled]), which makes the flux level susceptible to 
the effects of seeing variations between the SN observations and those of the flux standard 
star; some observations were also not obtained under photometric conditions. To produce 
an approximate absolute flux calibration, then, we computed synthetic photometry on our 
spectra and compared it with estimates of the V— or i?— band (depending on the spectral 
range of the spectrum) magnitudes of the SNe at the spectral epochs. Since none of 
our spectra cover an entire passband, we extended each spectrum's range by joining it with 
SN la spectra of similar age from the SUSPECT0 database^ scaled to match the flux of our 
spectra in regions of overlap. We then compared our synthetic photometry with estimates 
of the actual brightnesses of the SNe. Since no published late-time (i.e., beyond t > 100 
days) photometry exists for either SN 2005am or SN 2005cf, but high-quality early-time 



photometry does (ILi et al.ll2006l : iPastorello et al.l 120071 ). we derived approximate V and R 



magnitudes of the SNe at our spectral epochs by flrst extrapolating the early-time light 



curve s to day 200 through comparison with the light curves of SN 2003d u (IStanishey et al. 



20071 ). and then applying the late-time Type la decay-rates reported by iLair et al. J200d : 



AV = 1.46 mag 100 day^\ AR = 1.54 mag 100 day^^). We then adjusted each of our 
spectra by multiplying it by the scale factor needed to bring its synthetic magnitude to the 
estimated actual apparent magnitude. The only spectrum for which we did not apply this 
procedure was the day 298 epoch of SN 2005am, since the night was photometric, the flux 
of the SN from one exposure to the next was very consistent (indicating non- varying seeing 
throughout the observation sequence), and the spectrum of a nearby star (HD 79289) taken 
immediately after the SN observations was found to have a synthetic V magnitude within 
16% of the value reported by SIMBAD@ the scale factor for this epoch was thus taken to be 
the value needed to place HD 79289 on a correct absolute scale. The scale factors applied 
to all spectra are given in Table [2l Given the uncertainties inherent in our technique, the 
absolute flux calibrations are probably accurate to ~ 25%, in general. 

We next removed the redshift and rebinned each rest-frame spectrum to 3 A bin~^, the 
approximate spectral resolution (see Tabled]). For SN 2005am in NGC 2811, we adopted 
the NASA/IPAC Extragalactic Database (NED) recession velocity of 2,368 km s~^, and 
for SN 2005cf in MCG-01-39-3 we adopted the NED recession velocity of 1, 937 km s-\ 



see 



http://bruford.nhn.ou.edu/ ~suspect/indexl.htinl 



The specific spe ctra used for this pu rpose include SN 1998bu (day 329: ICappellaro et al.l 120011 ). 



SN 2003cg (day 385: lEhas-Rosa et al.lliooih . and SN 2002bo (day 375; 



Benetti et al 



20041) . 



see |http: / / simbad.u-strasbg.fr/simbad/ 
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We then searched for narrow Ha emission within ±1000 km s"^ (±22 A) of Aq = 6563 A 
in each rescaled, rest-frame spectrum by smoothing t he spectrum with a second-order 
Savitsky-Golay smoothing polynomial ( Press et al.lll992 ) of width ~ 100 A, differencing the 
smoothed and unsmoothed spectra, and then examining the residuals for narrow emission 
near Ha. This procedure yielded no evidence for narrow Ha emission in any spectrum. 
Indeed, no unexpected narrow features at any wavelength (including spectral regions near 
[O I] AA6300, 6364 and [Ca II] AA7291, 7324, which might be expected to produce features 
similar to Ha at these epochs from stripped solar-abundance material) were found. 

To determine the greatest strength (i.e., equivalent width) of a feature that could have 
remained undetected a t the l ocation of Ha in each spectrum, we adopted the procedure of 
Leonard &: Filippenkd (120011 ) . who report the 3o" lower bound on the equivalent width of 
an undetected feature in a spectrum to be: 



1^a(3(t) = 3AA M^Wune/^\ VVB, 



(1) 



where AA is the width of a resolution element (in A), A J is the la root-mean-square 
fluctuation of the flux around a normalized continuum level, is the full-width at half- 

maximum of the expected line feature (in A) and B is the number of bins per resolution 
element in the spectrum. The values for AA, A/, and B, along with the computed values 
of Wa(3o"), are reported in Table [2J The width of the line feature, Wune, was assumed to 
be 22 A. 

While m any studies have pre dicted the existence of narrow Ha in nebular spectra of 
SNe la, only iMattila et al.l (l2005l ) gives a quantitative estimate of its expected strength. 
Thus, to translate our detection thresholds into estimates of the maximum amount of 
solar-abundance material that could have given rise to an Ha line too weak to have been 
detected, we rely on this study alone. 



In the iMattila et al.l (120051 ) study, the code described by iKozma et al.l (120051 ) is modi- 
fled to artiflcially include varying amounts of solar abundance mater ial in the inner (±1, 00 km 



s ^) region of the ejecta. The resulting spectra were computed by IMattila et al.l (120051 ) in 



a time-dependent, one-dimensional, spherically symr aetricu manner for an SN la 380 days 
after maximum light. From the parameters given by IMattila et al.l (see their Figure 6 and 
associated discussion), we estimate that 0.05 Mq of solar-abundance material produces a 
luminosity in the peak of the Ha line of 



3.36 X 10^^ erg s"^ A" 



Translating predicted Ha luminosity into observed Ha strength requires estimates of 



^Note that while the expected asymmetr y in the location of th e stripped material would undoubtedly 



alter the shape of the resulting line profile, Mattila et al. ( 20051 ) conclude that it should not affect the 



total strength of the line, and so the assumption of spherical symmetry is deemed to be acceptable for the 
purposes of their study (and, hence, ours). 
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the distance and extinction of the SNe. For SN 2005am, we adopt the NED Hubble flow 
{Hq = 73 km s ~ ^ Mpc~^) distance o f 36.8 ± 2.6 Mpc based on the 3k CMB rest frame 



as determined 
determined by 



Li et al 



J2006h. 



Fixsen et al.l (119961 ). and the total extinction of Ay = 0.27 ± 0.08 mag 



For SN 2005cf, we use the distance and extinction values 
adopted bv iPastorello et all tOO^ oi D = 31.77 ±4.8 Mpc and Ay = 0.32 ±0.03 mag. For 
0.05 Mq of solar- abundance material 380 days after explosion, we derive expected Ha line 
fluxes at the profile's peak of 1.72 x 10"^^ and 2.23 x 10"^^ erg s"^ A"^ cm^^ for SN 2005am 
and SN 2005cf, respectively. 

To derive the equivalent width that such features would have had in our observed 
spectra, we approximated the Ha emission as a Gaussian with a la width of 11 A, and 
placed the resulting profile on top of the estimated "continuum" level at A = 6563 A (see 
Fig. [T]). We then measured the equivalent width of these theoretical profiles within the 
region ±11 A from Ao = 6563 A [the same spectral range used to derive Wx{3a) earlier] 
in each spectrum. The calculated values for W\{0.05Mq) are given in Table p Si nce the 
Ha emission is clearly a time-dependent phenomenon and the iMattila et al.l (120051 ) model 
is specifically for day 380, our estimate of the strength of the line in the day 298 and day 
267 spectra of SN 2005am and SN 2005cf, respectively, are quite crude. Fortunately, since 
the Ha emission is powered primarily by gamma-ray deposition (IMattila et al.ll2005l ). and 
the optical depth to gamma-rays in the central, hydrogen rich region should be higher at 
earlier times, we would expect the Ha line to be even stronger at earlier times than it is 
on day 380 (assuming, of course, that the outer, iron-rich ejecta are transparent at earlier 
times, and allow us to see it; see § [3]). Thus, our calculated value of Wa(0.05 Mq) for 
epochs earlier than 380 days should serve as a conservative lower limit for these spectra. 

Finally, we assume a linear scaling between the amount of stripped material and the 
expected equivalent width of the Ha emission line (an approximation justified through 
examination of Figure 6 of IMattila et al.l |2005| . where theoretical profiles are shown for 
0.01 and 0.05 Mq) to arrive at our final estimate for the upper limit on the amount of 
solar abundance material that could have remained undetected at each spectral epoch. We 
calculated this limit according to: 



M{Mq) < 



Wx{3a) 



Wa(0.05Mo) 



X 0.05. 



(2) 



Final results are listed in Table [21 and can be summarized succinctly: We restrict the 
amount of solar-abundance material that could have evaded detected to be < 0.01 Mq 
at all spectral epochs except day 384 of SN 2005cf, for which the formal upper limit is 

0.02 Mq. 
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Discussion 



When c oupled with the mode ling; r esults of iMattila et al.l (120051 ) and the mass-stripping 
estimates of [Marietta et al.l (120001 ) and lMeng et al. J2007I ). our restrictions on Ha emission 
place strong constraints on the progenitor systems that could have given rise to SN 2005am 
and SN 2005cf. In short, they rule out all hydrogen- donating companions close enough 
to the WD to have been experiencing RLOF at the time of explosion (i.e., a/R < 3; 
see § [1]). Under the favored single-degenerate scenario, this leaves only widely separated 
systems in which the companion (i.e., a red giant) donates matter through a strong stellar 
wind, as a viable option. For the symbiotic case, if we assume that the same scahng 
between the am ount of stripped material and the orbital separation holds for red giants 
as was found by [Marietta et al.l (120001 ) for main-sequence sec ondaries (see § [T| note that 
red giant companions were not analyzed by lMeng et al.l 120071 ). then we conclude that the 
mass-donating red giant must be > 10 A.U. from the WD to be consistent with our Ha 
nondetections. For a 1.4 Mq primary and a 1.0 Mq secondary, this corresponds to an 
orbital period of > 20 yr. 

Our main conclusion — that the most likely progenitors of SN 2005am and SN 2005cf 
were w idely separated symbiotics — stands in evident contrast to the work of lPanagia et al. 



(120061 ) on t he search for, and su bsequent nondetection of, radio emission from SNe la. In 
their study, iPanagia et al.l (120061 ) set upper limits on mass-loss rates of ~ 10~^ Mq yr~^ for 
the progenitor systems. Since thi s is of the same order of magnitud e as the observed mass- 



loss rates from symbiotics (e.g., lJung &: Led 12004 : ICrowleyl l2006l ). iPanagia et al.l (120061 ) 



specificall y rule out synabiotics as potential progenitors of SNe la. However, the inferences 
drawn by iPanagia et al.l (120061 ) rely on the assumption that SN la radio light curves behave 
identically to those of SNe Ib/c. Such an assumption is necessary because the correlation 
between wind density and radio luminosity for SNe la cannot be calculated from theory at 
this point. In the absence of theoretical or observationa l supp ort for such an assumption 
the very l ow mass-loss rates claimed by iPanagia et al.l (120061 ) cannot yet be considered 
definitive (IHughes et al.li2007i ). The radio non-detections may therefore still permit the 
symbiotic scenario. Along th ese same lines, the investigation of X-ray non-detections of 
SNe la by iHughes et al.l (120071 ) provides mass-loss limits of only ~ 10~^ Mq yr~^, which are 
not stringent enough to rule out symbiotics. While a low-density circumstellar environment 
is clearly favored by both the radio and X-ray non-detections of SNe la, it seems premature 
to rule out progenitor classes on this basis at the present time. 

There are additional caveats to our conclusion that must be given. Of primary signifi- 
cance is the fact that its robustness depends critically on the validity of the theoretical mod- 
eling of just a few studies. Taken together, these st udies predict that a substantial amount 
of stripped hydrogen should exist at low ve locities (iMarietta et al.ii2000l : iMeng et al.li2007i ) 



and emit a strong Ha feature at late times ( Mattila et al. [2OO5I ). There are, however, alter- 



native explanations for the observed lack of hydrogen that must be thoroughly investigated 
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to build confidence in the inferences tliat can be drawn from our study. We now critically 
evaluate four such alternatives, and point out areas ripe for additional theoretical study. 

(1) The hydrogen is "hidden" behind an opaque screen of faster-moving ejecta. Posed 
simply: Are the fast-moving, iron-rich ejecta truly optically thin at late times? Since per- 
mitted Fe II lines dominate the underlying spectrum of SNe la near 6563 A, significant 
opacity in these lines co uld provide a n "iron curtain" that shields the hydrogen emission 
from our view. Indeed, iKozma et al.l (120051 ) state that several Fe II lines do remain op- 
tically thick in thei r models even at late times. On the other hand, the simulations of 



Matti 



els of 



Kozma et al 



a et al.l (120051 ) , whi ch predict the strong Ha feature, are exphcitly based on the mod- 



(120051 ) . and so it seems reas onable to infer that the specific Fe II lines 
near Ha do not provide significant opacity in the lKozma et al.l (|2005l ) models. Independent 
confirmation of the transparency of the outer ejecta of SNe la at late times is needed. 

(2) There is far less hydrogen entrained in the ejecta than current models predict. By 
predicting up to a factor of four reduction in the amount of stripped mass compared with 
the models of iMarietta et al.l (120001 ). iMeng et al.l (120071 ) demonstrate the potentially large 
effect that binary evolution can have on the amount of stripped material. Might other 
aspects of binary evolution be at work t o further re d uce t he stripped mass? A possibly 
important mechanism not considered by iMeng et al.l (120071 ) is the outcome that a strong 
"accretion wind" , blown from the WD during the mass-accretion ph ase, could have on the 
secondary's envelope. Accretion winds were originally conceived by iHachisu et al.l (119961 ) 
as a means to stabilize the mass-transfer process and al l ow th e sys tem to avoid a common 
envelope phase. Further investigation by lHachisu et al.l (119991 ) and lHachisu fc Katol (120031 ) 
revealed that the winds should also heat and ablate a si gnificant amount of m aterial from a 
close secondary's envelope. From these considerations, iHachisu et al.l (119991 ) estimate that 
the wind could strip mass from such a companion at rates as high as ~ 10~^ Mq yr~^ for 
a period of up to ~ 10^ yr prior to the explosion. If such a process actually precedes an 
SN la explosion, it would clearly leave the secondary's envelope in a vastly different state 
from the models considered to date. In fact, such secondaries could potentially lose their 
entire hydrogen envelope and evolve to become he lium-donators, whi ch would naturally 
account for a lack of hydrogen in late-time spectra (jSranch et al.lll995l ). 



There are both theoretical and observational objections to strong WD winds playing a 
large roll in stripping ma ss from the s econdary star in SN la progenitor systems, however. 
On the theo retical front. iLivid (120011) questions the high efficiency of the mass stripping 
proposed by IHachisu et al.l (119991 ) . and also points out that at high accretion rates, the 
bulk of the WD's wind may be strongly collimated in a direction perpendicular to the 
accretion disk rather than in the direction of the companion star. Such a scenario would 
result in little stripping. Obs ervational doubt on a ccretion winds playing a major role in 
pre-SN la evolution is cast by lBadenes et al.l (120071 ). who find that such an optically thick 
outflow from a WD's surface should excavate a large, low-density cavity around the system. 



- 10 - 



Such a large evacuated region is incompatible with the known observational properties of 
SN la remnants in the Galaxy, the Large Magellanic C loud, and M31. It thus appears that, 
while such accretion winds may occur in nature (e.g.. iHachisu &: Katdl2003l ). they do not 



evidently precede a majority of SN la explosions. Further consideration of this and other 
potential mass-stripping mechanisms is certainly warranted. 



(3) The hydrogen gas is insufficiently powered to produce Ha emission. iMattila et al. 



(120051 ) conclude that the hydrogen-rich, central region of the young supernova remnant 
presents high optical depth to gamma-ray photons that originate from the radioactive decay 
of isotopes synthesized in the explosion. The gamma-ray trapping is then responsible for 
powering the Ha line emission. This mechanism is quite different from the situation in the 
outer, lower density regions of th e ejecta, whose optical eniission is predominantly powered 



by local positron deposition (e.g. JStritzinger fc Sollermanll2007l ). Since the production of a 



strong Ha line from even small amounts of solar-abundance material relies fundamentally 
on the inner ejecta's ability to trap gamma rays, the high opacity of the material, and the 
existence of sufficient numbers of gamma-rays to power the emission, should be confirmed 
by independent modeling. 

(4) The single- degenerate scenario is not responsible for these SNe la. A double- 
degenerate progenitor system, in which two WDs in a binary system coal esce due to the 
emission of gravitatio nal radiation and ultimately explode as an SN la (IWebbinkI 1 1984 : 
Iben fc Tutukovlll984j ). would naturally account for the absence of hydrogen in the spec- 
tra. Indeed, detection of any hydrogen would deal a death blow to the double- degenerate 



scenario. Difficulties with the double-degenerate scenario are well known (see iLivid 12001 
for a thorough review), however, and include both observational (e.g., of the ~ 120 double- 
degenerate systems kn own, none have a total mass in excess of the Chandrasekhar limit; 
Napiwotzki et al.l 120041 ) and theoretical (e.g., the merger process seems m ore likely to lead 
to co llapse to a neutron star than to thermonuclear explosion as an SN la; ISaio fc Nomoto 
19981. and references therein) objection s. Nonetheless, theoretical loopholes remain (e.g., 
Piersanti et al.l l2003l : lYoon et al.l 120071). and the observati onal sample of nearby double- 
degenerate systems is not complete (INelemans et al.ll2005l ). Since the double-degenerate 
channel is a potential "silver-bullet" that can explain both the lack of hydrogen in SNe la 
as well as the occ urrence of SNe la in both old and young star-forming systems (e.g.. 
Branch et al.l 119951 ). it requires continued theoretical and observational attention. 



4. Conclusions 

We obtained five deep, moderate-resolution, nebular-phase spectra of two SNe la 
(SN 2005am and SN 2005cf) in order to search for narrow Ha emission that would betray 
the existence of material stripped from the envelope of a mass-donating stellar companion 
to the exploding WD. No such emission is detected in either object at any epoch. From the 
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models of iMattila et al.l (120051 ) . we establish upper limits of 0.01 Mq of solar abundance 
material in the inner ejecta of both objects, which are the tightest constraints yet estab- 
lished by such s t udies . Our non-detections of Ha, coupled with the mass-stripping results of 
Marietta et al.l (120001 ) and iMeng et al.l (120071 ). rule out all hydrogen- donating companions 
close enough to the WD to have been experiencing RLOF at the time of explosion for these 
events. Additional theoretical work is needed in several areas to buttress this conclusion, 
including most critically verification of the transparency of the outer, more rapidly mov- 
ing ejecta that could potentially block Ha photons from escaping from the inner region. 
Bearing this caveat in mind, we propose that symbiotics are, at this time, the most likely 
progenitor class that remains consistent with these data. 

Definitive proof of the identity of the progenitor system(s) that gives rise to SNe la 
remains elusive, and it must be admitted that our conclusion, which is based on the lack 
of a detection, is not as satisfying as one based on a detection. Should future modeling 
efforts prove unable to "hide the hydrogen" for even widely separated binaries, then the 
continued viability of the single-degenerate, hydrogen- donating progenitor will require that 
Ha, no matter how weak, must ultimately be detected. 
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Fig. 1. — Late-time spectra of two SNe la, with day since B maximum indicated. The 
spectra are displayed at 3 A bin~^, the approximate resolution at Ha. The expected 
strengths of the Ha lin e resulting fro r n 0.05 Mq of solar- abundance material according to 
the day 380 models of iMattila et al.l (120051 ) are shown as dot-dashed lines in the insets; 
note that since the Ha emission is a time-dependent phenomenon, the estimated strength 
of these lines in the day 298 and day 267 spectra of SN 2005am and SN 2005cf, respectively, 
are only approximate. 



Table 1. Journal of Spectroscopic Observations 









HJD 




Range^ 


Resolution^ 


P.A.= 


Par. P.A.'' 




Flux 


Seeing^ 


Slit 


Exposure^ 


UT Date 


Object 


Day" 


-2,400,000 


Telescope^ 


(A) 


(A) 


(dcg) 


(dcg) 


Air MassS 


Standard^^ 


(arcscc) 


(arcscc) 


(s) 


2005 Doc 31.52 


SN 2005am 


298.02 


53736.02 


KI 


5860-7160 


3.0 


150 


149-4 


1.24-1.32 


HD84 


1.3 


1.0 


7,200 


2006 Mar 6.59 


SN 2005cf 


267.09 


53801.09 


GN 


5520-7628 


3.1 


55 


135-22 


1.12-1.33 


HZ44 


0.7 


1.0 


10,800 


2006 Mar 24.15 


SN 2005am 


380.65 


53818.65 


GS 


5150-6920 


3.2 





117-156 


1.04-1.34 


LTT1020 


0.7 


1.0 


10,410 


2006 Jun 1.46 


SN 2005cf 


353.96 


53887.96 


KI 


5750-7430 


2.4 


55 


21-58 


1.15-1.70 


BD26 


1.1 


1.0 


11,100 


2006 Jul 1.39 


SN 2005cf 


383.89 


53917.89 


KI 


5796-7496 


3.4 


55 


47-61 


1.32-1.93 


BD26 


1.0 


1.5 


6,600 



Note. — All Keck observations were obtained with the 900/5500 [number of lines mm "'"/blaze wavelength (A)] grating + GG495 order blocking filter; all Gemini observations were carried out with 
the R831_G5322 grating + OG515_GO330 order blocking filter. Gemini observations were conducted under programs GS-2006A-Q-28 {SN 2005am; PI: Leonard) and GN-2006A-Q-27 (SN 2005cf; PI: 
Leonard) . 

^D ays since estimated dates of maximum B brightness. SN 2005am: 2005 March 8.5±1, HJD 2, 453, 438±1 jBrown et al|2005l ). SN 2005cf: 2005 June 12.5it0.3, HJD 2, 453, 534.0±0.3 l lPastorello et al.| 
[2003). 

^K I ^ Keck I 10 m/Low-Rcsolution Imaging Spectrometer + polarimcter {LRISp; lOke et al1ll995l ): GN(S) ^ Gemini North (South) 8 m/Gcmini Multi-Object Spectrograph (GMOS; iHook et all 
f200i) . 

'^Wavelength range of the calibrated flux spectrum. 

'^Approximate spectral resolution derived from night-sky lines near the nominal location of Ha in the observed spectrum. 
^Position angle of the spectrograph slit. 

Parallactic ang le iFiliDDenkdll982ll range calculated at the midpoint of each separate observation for each set of observations. 
^Airmass range calculated at the midpoint of each separate observation for each set of observations. 

^HD84 = HD 84937, BD26 ^ BD-|-26°2606 l lOke 8z GtI^|1983D : HZ44 ^ Hz 44 jMassey &z GrQnwall||199GD : LTT1020 ^ LTT 1020 jBaldwin &z StondfigS^ ) . Standard stars for observations made at 
the Keck Observatory were observed on the same night as the supernova observations; for the Gemini Observatory observations, the standard stars were observed on 2006 February 15 (HZ44, for the 
day 267 observation of SN 2005cf) and 2005 September 5 (LTT 1020, for the day 381 observations of SN 2005am). The observations of LTT 1020 were made as part of program GS-2005B-C-4 (PI: Dan 
Christlein) . 

^Average value of the full width at half maximum of the spatial profile for each set of observations, rounded to the nearest 0. 1. 

■I Combined exposure duration of all observations, in seconds. Observations at the Keck Observatory consisted of four separate, successive exposures, while observations at the Gemini Observatories 
consisted of six separate, successive exposures. 
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Table 2. Measured Values 







Scale 








Wa(O.O5M0)^ 




Object 


Day 


Factor*^ 




Ar 


(A) 


(A) 




SN 2005am 


298.02 


1.16 


1.00 


0.090 


2.19 


10.76 


0.01 


SN 2005cf 


267.09 


1.45 


0.97 


0.015 


0.38 


2.11 


0.01 


SN 2005am 


380.65 


0.45 


0.94 


0.144 


3.73 


21.51 


0.01 


SN 2005cf 


353.96 


1.22 


1.25 


0.055 


1.07 


8.82 


0.01 


SN 2005cf 


383.89 


1.02 


0.88 


0.164 


4.51 


13.10 


0.02 



Note. — All measurements made on rest-frame spectra rebinned to 3 A bin ^. 

'^Factor by which the original, reduced flux spectrum was multiplied to place it on an (approx- 
imate) absolute flux scale; see text for details. 

'^Number of 3 A bins per resolution element. 

■^IcT root-mean-square fluctuation of the flux around a normalized continuum. 

^Equivalent width of the strongest potential feature near Aq = 6563 A that could have remained 
undetected in the spectrum, derived using Equation [TJ 

'^Expected equivalent width of an Ha emission l i ne res ulting from 0.05 Mq of solar abundance 
material according to the models of lMattila et al.l (120051 ): see text for details. 

^Derived upper limit (rounded up to the nearest 0.01 Mq) on the amount of solar abundance 
material that could have remained undetected at each spectral epoch, derived according to Equa- 
tion [2l 



